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ABSTRACT 

With the recent interest in the measured and standard model values of the muon 
anomalous magnetic moment, a M , some confusion has arisen concerning our knowledge 
of the hadronic contribution to a^. In the dispersion integral approach to hadronic 
vacuum polarization effects, low energy contributions must be evaluated from data or 
in a model-dependent approach tested by data. At higher energies perturbative QCD 
has been used, sometimes in conjunction with data. The history of such evaluations 
is reviewed, and the prospects for further improvement are discussed. We conclude 
that not all published evaluations are on an equal footing or up-to-date. One must 
critically examine which, and how much information went into each analysis in order 
to determine which are more complete, and reliable. 



1 Introduction 

With the new result from the Muon (g — 2) Collaboration |], comparison with the 
standard model prediction for has received renewed scrutiny. This result represents 
the third measurement by the E821 collaboration 0, |J, and the first to approach the 
part per million (ppm) level of precision. That level of precision permits a stringent 
new test of the standard model, and a search for physics beyond it. However, before 
conclusions can be drawn, it is necessary to have a reliable standard model calcula- 
tion with sub-ppm accuracy for comparison with the experimental number. In this 
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paper, the status of QED and weak contributions are briefly reviewed and then we 
focus on the hadronic contributions. After an overview of all published evaluations 
of the hadronic contribution since 1985, we discuss several in detail, with an eye to- 
wards assessing both reliability, and whether the calculation is current or has been 
superseded as a point for comparison with experiment by more up-to-date studies. 
The standard model prediction for a M = (g^ — 2)/2 consists of three parts, 

a M (theory) = a^(QED) + a /t (Hadronic) + a^Weak). (1) 

Taking the value of a from the electron (g — 2) |4]], yields the total QED 
contribution || U 

a M (QED) = 116 584 705.7(2.9) x KT 11 , (2) 

which is dominated by the first-order (Schwinger) term a/2ir but is calculated (or 
estimated) through 0(a/n) 5 . The uncertainty is very small and should not play a 
role in comparisons with experiment. 

The weak contribution through second order is @, §, [|, |i~0H 

a^(weak) = 152(4) x KT 11 (3) 

contributing about 1.3 ppm of a M (assuming a 150 GeV Higgs mass). The 3-loop 
electroweak leading logs have been estimated to contribute +0.5 x 10 -11 . That small 
effect is safely covered by the uncertainty in eq. @. 

Although QED and electroweak effects now appear to be well under control, there 
have been some sizeable shifts in their predicted values over the years due to error 
corrections and improved higher order calculations. In Table p], we illustrate changes 
in the theoretical prediction for a M that have occurred since a summary was given in 
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While the QED and weak contributions are well described perturbatively, the 
hadronic contribution cannot be completely calculated from perturbative QCD, but 
must instead be determined in part by using data from the cross-section 



for [15, 16, 17, 18, 19, 20, 21 



e + e — > Hadrons (4) 

in conjunction with a dispersion integral. The uncertainty in those data, particularly 
at low energies, largely determines the error in the Standard Model prediction for 
a^. More recently, data from hadronic r decays have also been used along with 
information from perturbative QCD at relatively low energies [f22|, ^3, |24| to reduce 
uncertainties. 
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Quantity 


1550 Value 
(xlO 11 ) 


2001 Value 
(xlO 11 ) 


Change (xlO 11 ) 


a™ d (vac. pol 1) 

a^ ad (vac. pol 2) 
a Had j light by light ) 

a™ (1 loop) 
af w (2 loop) 


116 584 695.5(5.4) 
7 068(59)(164) 
-90(5) 
49(5) 
195(10) 


116 584 705.7(2.9) 
6 924(62) 
-100(6) 
-85(25) 

195 
-43(4) 


+10.2 
-144 
-10 
-134 


-43 


af 1 (total) 


116 591 918(176) 


116 591 597(67) 


-321 



Table 1: Improvements in the theoretical calculation of a M from 1990 [[14]] to 2001. The 
major shifts were primarily due to errors in the earlier calculations, new calculations 
of higher order effects, improved e + e~ — > hadrons and tau data, and additional 
utilization of perturbative QCD. 



2 The Hadronic Contributions 

In this section we describe the basic issues and summarize the various evaluations of 
the hadronic contributions to the muon anomaly, including the higher-order 3-loop 
effects. Then we discuss in some detail one of the data-driven analyses in order to 
illustrate the relative importance of the various energy regions in the evaluation of the 
leading hadronic contribution, a M (Had; 1), and the main sources of the errors on it. 
Finally, we discuss some specific details of several of the calculations. In particular, 



it is hard to ignore the recent assertion |25j that we currently have only modest 
knowledge of a M , and all calculations are of equal merit. In this latter section we will 
refute the claims made in that paper. 



2.1 The Leading Hadronic Contribution 

The leading hadronic contribution comes from the vacuum polarization diagram 
shown in Fig. ffl(a). Because the loop integration involves low energy scales near the 
muon mass, the contributions cannot be calculated from perturbative QCD alone. 
At higher loop momenta perturbative QCD becomes applicable, and it is common in 
the evaluations of a M (Had; 1) to effectively switch from data to QCD at some energy 
scale. 

Since the leading order hadronic vacuum polarization contribution is derived pri- 
marily from data, its value continues to be improved by new e + e~ cross-section mea- 
surements. In that regard, the CMD2 collaboration at the VEPP2M collider in 
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Figure 1: (a) The leading hadronic vacuum polarization contribution, (b-d) Examples 
of higher order (a/ir) 3 contributions derived from the hadronic vacuum polarization. 



Novosibirsk has collected substantial data from the 2tt threshold up to y/s = 1.4 
GeV, and is preparing a publication on these data [[Ti l. Similarly, the BES collab- 
oration in Beijing has recently measured R(s) from yfs = 2.0 GeV up to 5 GeV 
[|j~2]j . Data from hadropro duct ion (see Fig. 0(a)) can be related through a dispersion 
relation to the first order hadronic vacuum polarization of Fig.|l](a). 

A second way in which the hadronic contribution can be improved is through 
hadronic r decays to vector final states. The relevant diagram is shown in Fig. 
|2|(b) where the weak charged current can be related to the isovector part of the 
electromagnetic current in Fig. 0(a) through the conserved vector current (CVC) 
hypothesis plus the additional requirements of isospin conservation and the absence 
of second-class currents (which is the case for the Standard Model). While the weak 
t decay proceeds through both vector and axial vector weak-currents, the final states 
with an even number of pions (even G-parity) are the ones relevant for (g — 2), since 
decays to these final states go exclusively through the vector current if there are no 
second class currents. Of course, isospin violating effects of order a few percent must 
be properly accounted for. 




(a) 

Figure 2: (a). The hadroproduction process which enters the dispersion relation, 
(b) Hadronic r decay. 
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The measurements of 



R 



^totle^e 
o"tot(e+e- 



are used as input to the dispersion relation, 
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Since there is a factor of s~ 2 in the dispersion relation, values of R(s) at low energies 
dominate the calculation of a^(Had; 1). 

Here, we should point out that R is often not directly measured experimentally. 
In those cases where the cross-section for e + e~ is determined using some other nor- 
malization, careful subtractions for initial state radiation, vacuum polarization etc. 
have to be applied to the data whereas most such effects would naturally cancel in 
the ratio R. In Table [| the published evaluations of a M (Had; 1) since 1985 are given. 

In carrying out the dispersion integrals, there are two ways of integrating the 
data. In B85 and KN085 the e + e~ data were fit to models such as a Gounaris- 
Sakurai resonance parameterization or extensions thereof, and these theoretical 
curves were then utilized to obtain R(s) for the dispersion integral. EJ95 [19| were 
the first to employ a model-independent trapezoidal integration of data. BW96 |pOfl , 
also used a trapezoidal integration and took into account the correlations between the 
systematic errors in the data. Below we employ one of these two model independent 
evaluations, which are in excellent agreement with one another, as a benchmark to 
discuss other evaluations of the hadronic contribution. 



2.2 Higher Order Hadronic Contributions 

The higher-order (3-loop) hadronic contributions come from the hadronic vacuum 
polarization diagrams in Fig.|l|(b-d) and the hadronic light-by-light scattering shown 
in FigfJ 

The higher order hadronic contributions in Fig|lj(b)-(d) were most recently eval- 
uated by Krause p9| . The set of diagrams represented by Fig. |](d) give a dominant 
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Ref. 



,(Had; 1] 



xlO 



a M (Had; 1) in ppm 



Comments 



B85 |g 
KN085 [0 
CLY85 HT 
MD90 1 



MD90 p 



EJ95 P] 



BW96 |0| 
AY95 fl2H 



ADH98 |2| 



6840 (110) 
7070 (180) 
7100 (115) 
7048(115) 
7052(76) 
7024 (153) 
7026 (160) 
7113 (103)* 
6950 (150) 



58.6 
60.6 
60.9 
60.5 
60.5 
60.3 
60.3 
60.0 
59.6 



0.9) 

1.54) 

0.9) 

1.0) 

0.65) 

1.4) 

1.4) 

0.9) 

1.29) 



primarily e + e~ data 
primarily e + e~ data 
QCD, theory and some e + e" 
e + e- and QCD 
e+e" and QCD 
primarily e + e~ data 
primarily e + e~ data 
QCD, theory and some e + e~ 
primarily e + e~ data 



ADH98 |2| 
DH98a @ 



DH98b 24 



7011 (94) 
6951 (75) 

6924 (62) 



60.1 
59.6 



0.8) 
0.6) 



59.4 (.53) 



e e~ + t data 

e + e _ , t and perturbative QCD 
at energies (E > 1.8 GeV) 
e + e~, r and perturbative QCD 
and QCD sum rule constraints 
at low energy 



Table 2: The first-order hadronic vacuum polarization contribution to (g — 2) obtained 
by a number of different authors. (Earlier evaluations with much larger uncertainties 
are not exhibited.) 

*The value of a M (Had; 1) given in the abstract of |2lH ([hep-ph/ 9509378]) does not agree 
with the value given in the conclusions section of text. Perhaps this confusion was 
corrected in the 1998 published version which is not available to us. We assume that 
the value quoted in [^j is the value to take from |?TJ. Similarly, in [17] a second 
method for evaluating the p contribution gave a somewhat smaller result 7045 which 
is not illustrated in the table. 



negative contribution. Collectively, Krause found —101(6) x 10 n . That value was 
slightly updated by Alemany, Davier and Hocker E3] to 



a^Had; 2) = -100(6) x 10" n (0.86 ± 0.05)ppm (9) 



which supersedes the earlier value |nj . The difference between the earlier evaluation 
and the more recent calculation is attributed to the use of the full kernel function in 



the new calculation p9[, and more up-to-date hadronic data. 
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u -11 
r -85 (32) X 10 



Figure 3: The hadronic light-by-light scattering contribution. 



2.2.1 Hadronic Light-by-Light Contribution 

The hadronic light-by- light contribution shown in Fig.|]was first calculated by KN085 
Unlike the hadronic contributions discussed above, this contribution cannot be 
determined from data, so one is dependent on a model calculation. Einhorn |3(| 



pointed out that the vector meson dominance model used in KN085 did not satisfy 
the Ward-Takahashi identities which are required by electromagnetic gauge symmetry. 
Furthermore, part of the calculation was in error. 

There have been three recent calculations of the hadronic light-by-light scattering 
contribution, one by Hayakawa, Kinoshita and Sanda (HKS) one by Bijnens, 
Pallante and Prades (BPP) ||32|| , and a follow-up improved calculation by Hayakawa 
and Kinoshita (HK) |33|. Along the way, low energy theorems were developed by 
HayakawaQ for the (p-wave) V° — n scattering amplitude (V° = p°, u or 0) [ |34[ . 

The situation seems to have converged as far as is possible without a first princi- 
ples or lattice calculation of the four-point function. Since the light-by-light contribu- 
tion currently can only be obtained from calculation, the give and take between two 
groups, and the mutual checking of the other's calculation was invaluable in resolving 
the magnitude and uncertainty on this correction. In the final analysis, it may be the 
uncertainty of this contribution which provides the ultimate limitation on the stan- 
dard model prediction for a Ai (Had). Indeed, subtle cancellations in the calculation 
deserve further study. Fortunately, the current level of uncertainty appears suitable 
for the final goals of E821. 



Following [22], we take the average of (BPP) and (HK), 



a M (Had; lol) = -85(25) x KT 11 (-0.72 ± 0.21) ppm (10) 



The total higher order hadronic correction from (|9]) and ( |10] ) is 
a M (Had; Higher order) = -185 ± 26 x 10~ n 



Ref. p4| the behavior in the limit of low pion momentum is studied, and it is shown that the 
s-wave component must vanish in the chiral limit under quite general assumptions. 
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3 The Hadronic Contribution Compared with a^ xp 
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Figure 4: A comparison of the calculated values of a A1 (Had; 1) along with the difference 
between the measured value of a M less the QED, weak and higher order hadronic 
corrections. 



The recent result reported by E821 is 

a e * p = (116 592 020 ± 160) x 10" 11 (11) 

where the systematic and statistical errors have been added in quadrature. When 
combined with the previous measurements, one gets 

< a M >cx P = (116 592 023 ± 151) x lO -11 (12) 

Since the QED and electroweak contributions are not in question, we can subtract 
those calculated values from the experimental number. The resulting quantity rep- 
resents the hadronic contribution plus any contribution from new physics. Next, we 
subtract off the higher order hadronic contributions given above and obtain 
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a M (Had; 1) + a M (New?) = 7350(153) x KT 11 . (13) 

This number is to be compared with the values of a M (Had; 1) illustrated above in 
Table A graphical comparison is given in Fig. 

This compilation represents a completely uncritical selection of published evalu- 
ations. Nevertheless, all evaluations agree with each other and deviate in varying 
degrees with the value suggested by experiment. We discuss some of these in detail 
below. Because of the common data sets which went into most of these evaluations, 
the errors on these evaluations of a M (Had; 1) are highly correlated. 



4 The Different Evaluations of a M (Had; 1) 

We now give a detailed critique of the evaluations of a^(Had; 1). Rather than starting 
the discussion with B85, KN085, or CLY85, all of which used a parameterization for 
the p resonance we first consider the model-independent analyses of EJ95 and BW96. 

While the analyses of e + e~ data by EJ95, BW96 and ADH98 obtain essentially 
identical final results, there are differences on several points between the three analy- 
ses. In EJ95 two results are quoted. These are given in Table 3a, and Table 3b of ref. 
IS]. It is the former which is directly comparable to the other analyses. The so-called 
"renormalization group improvement" result quoted in EJ95 (Table 3b) represents an 
effort by those authors to incorporate the higher order contribution of Fig. |l|(b) into 
the lowest order hadronic vacuum polarization. Unfortunately, that approach is in- 
valid and should be discounted. The fact that it is quoted in the abstract of EJ has 



caused some recent confusion [25]. The authors (EJ) themselves have not continued to 
advocate the Renormalization Group improved approach, but instead have accepted 
the conventional method of computing higher order hadronic corrections separately 



as done by Krause.[|29] 



Although both EJ95 and BW96 do a model independent analysis, there is a dif- 
ference in how errors and experimental data are combined. Up to 2 GeV, the data 
have traditionally been published as exclusive cross sections. Above 2 GeV they are 
published as the inclusive cross section ratio R(s). It is necessary to integrate over 
energy, sum over modes (below 2 GeV), and combine results from different detectors. 
EJ95 computed an error weighted average of -R(s) over all experiments, but in so 
doing, they lost the information on correlations between the errors over energy. 

BW96 combined the data in a way which permitted the correlations to be included, 
and they also added a scale factor when combining results from different experiments, 
thereby handling experiments which do not agree in a conservative manner. 

While the EJ95 method gives up the information on the correlations over energy, 
and BW96 explicitly includes them, the final results are so close that there is no prac- 
tical difference between the two methods or results. We conclude that the correlations 
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are not important. Furthermore, when the better data which are now available from 
BES and other e + e~ facilities are included in the future, the issue becomes irrelevant. 

4.1 The Contributions to a M (Had; 1) from Different Energy 
Regions 

It is important to understand the source of the errors on a M (Had; 1) how they have 
recently been improved and will be further improved over time as new data become 
available. In Table [| the contributions to a A1 (Had; 1) from different energy regions 
in the dispersion integral found in 1996 by Brown and Worstell are listed with their 
uncertainties. 



Energy Region 


a^Had; 1) 


%of 


<5tot 


"frac 


(GeV) 


(xlO 11 ) 


a^Had; 1) 


(xlO 11 ) 


(ppm) 


<r(e + e — > Hadrons) y/s < 1.4 


6113.32 


87 


149.97 


1.29 


a(e + e~ -> Hadrons) 1.4 < ^/s < 2.0 


324.66 


4.6 


24.96 


0.21 


R{s) a 2.0 < y/s < 3.1 


283.74 


4.0 


35.51 


0.30 


R(s) b 2.0 < v^< 2.6 


204.8 


2.9 


31.47 


0.27 


R(s) b 2.6 < v/i< 3.1 


78.93 


1.1 


13.99 


0.12 


J/ifj (6 states) 


90.47 


1.3 


9.69 


0.08 


T (6 states) 


1.09 




0.13 




QCD 3.1 < v/i< oo 


213.01 


3.0 


3.71 


0.03 


Subtotal v/i < 3.1 + J ftp + T (no QCD) 


6813.28 


97 


160.22 


1.37 


Subtotal y/s < 1.4 


6113.32 


87 


149.97 


1.29 


Subtotal 1.4 < v/i < 3.1 


608.40 


8.7 


55.51 


0.48 


Subtotal t/s > 1.4 (Includes QCD) 


912.97 


13 


56.48 


0.48 


Total hadl 


7026.29 


100 


160.25 


1.37 



Table 3: The 1996 contributions to a^(Had; 1) from the various energy regions 
with their total errors. The systematic errors are twice the statistical errors in the 
region from threshold to 2.0 GeV and from 2.6 to 3.1 GeV. From 2.0 to 2.6 GeV the 
systematic and statistical errors are about equal. (This table is based on Table IX. 



from Ref. pi). 



The largest contribution comes from threshold to 1.4 GeV. However, in the context 
of this model-independent analysis, one sees that if the error from the 2ir threshold 
to 1.4 GeV were eliminated completely, one would still have been left with an error 
of 8 t ot = 56.5 x 10 -11 , or 0.48 ppm. The recent Aleph r-decay data have been used 
to significantly improve the region below 1.4 GeV, and inclusion of the full LEP and 
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CLEO r data samples along with the new e + e~ Novosibirsk data will improve things 
further. 

The region between 2.0 and 2.6 GeV had been particularly problematic. In 1996 
there were only two experiments, BCF and 772, which have comparable total errors, 
but 772 has a large systematic error. This small region alone contributed an un- 
certainty of ~ 0.27 ppm to a^(Had; 1), while the entire region from 2.0 to 3.1 GeV 
contributed ~ 0.31 ppm error. (Without correlations these errors add in quadrature.) 
With the new data from BES [fT2]] , which agree with the QCD evaluation by Davier 
and Hocker, it appears that the problems in this energy region are solved and the 
uncertainty is significantly reduced. Hence, more current studies of a /J (Had; 1) have 
justifiably smaller uncertainties than earlier efforts. 



4.2 Brief Overview of the Evaluations of a /x (Had; 1) 

Before giving a more detailed discussion of the evaluations listed in Table |2], we pro- 
vide a brief perspective. Our point of view is that the knowledge of the hadronic 
contribution to is an evolving topic, that earlier analyses which represented the 
state of the art at one time, become outdated when new data and improved evalua- 
tions become available. 

When B85 and KN085 made their analyses, there were many unpublished data 
from Orsay which were not included. By 1995 these data were mostly published, 
so we believe EJ95, BW96 and ADH98 supersede the earlier studies. We defer the 
discussion of CLY85 and AY95 to the next section, but note in passing that one of 
the goals of MD90 was to improve the estimate of the theoretical errors presented in 
CLY85. 

The values of a^(Had; 1) obtained by EJ95, BW96 and ADH98 are quite con- 
sistent. It is interesting to note the large improvement obtained by ADH98 by the 
inclusion of the r-decay data. From Fig. ^ and Table ^| one can see that the addi- 
tion of the r-decay data raised the value of a /J (Had; 1) slightly, but to a value quite 
consistent with the earlier analyses. The use of QCD for ^/s > 1.8 GeV by DH98a, 
along with e + e~ and r-decay data lowered the value of a M (Had; 1) and significanly 
reduced the error, since the QCD prediction was systematically below the existing 
data (discussed above) in the energy region 2-3 GeV. This QCD prediction, which was 
done in advance of the recently reported R(s) measurement at Beijing W2 



is m ex- 



cellent agreement with the new data, which do not seem to suffer from the systematic 
problems of the older BCF and 772 data |13| . 



This excellent agreement between the QCD calculation and the new R(s) data 
gives one confidence in the evaluation of a M (Had; 1) presented in DH98a. Since the 
input to this evaluation contains much more data than the earlier evaluations, and 
the theory input seems to be justified by the new Beijing data, we believe that one 
must at least take DH98a as the best evaluation of a M (Had; 1) up to that point. In 
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DH98b, the same authors use QCD sum rule constraints at lower energies to further 
improve the uncertainty on their evaluation of a At (Had; 1), and it is this last value 
which is used in Ref. |l[ to compare the experimental and theoretical evaluations. 
That final improvement is perhaps a more controversial improvement. However, we 
note that DH98a and DH98b do not significantly differ. 

To illustrate the degree of improvement in a^(Had; 1), we give in Table 4 contribu- 
tions from different energy regions as obtained by Davier and H6cker.|24| Although 



(GeV) 


a M (Had; 1) x 10 11 


2m« - 1.8 


6343 ± 60 


1.8-3.7 


338.7 ±4.6 


3.7-5 + ip(ls, 2s) 


143.1 ± 5.4 


5-9.3 


68.7 ± 1.1 


9.3-12 


12.1 ±0.5 


12-oo 


18.0 ±0.1 


Total 


6294 ± 62 



Table 4: Contributions to a^(Had; 1) from different energy regions as found by Davier 



and H6cker.[24 



the energy divisions in Tables 3 and 4 do not coincide, one can see indications of 
improvement throughtout Table 4. Particularly significant is the error reductions for 
2m n < y/s < 1.8 GeV due to tau data as well as for y/s > 1.8 GeV due to the use of 
perturbative QCD (now confirmed by BES data). 



4.3 Detailed Discussion of CLY85 and AY 

We now discuss some of the specific issues raised by Ref. [p5fl . Before undertaking 
this discussion, we note that the evaluations of a M (Had; 1) in CLY85 and AY95 are 
somewhat high, but consistent with other analyses, and for some of the same reasons 
that B85, KN085, EJ95, etc. are now considered to be obsolete, these two should also 
be so considered. We also recognize that CLY85 represented the first attempt to use 
perturbative QCD down to low energies in the evaluation of a /i (Had; 1). While it took 
some time to be accepted, the use of QCD, along with r-decay data has resulted in 
substantial improvements [^, ^4|. Nevertheless, we feel compelled by the confusion 



which has been recently generated by Ref. |25|] to detail why we feel that CLY85 and 



AY95 should not be used in comparison with experiment, at least not with the same 
conviction as later improved determinations. In fact, we argue that both should be 
updated before a serious comparison with more up-to-date approaches can be made. 
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Our first criticism of CLY85 and AY95 involves their treatment of the p resonance. 
Both employ a parametrization of the pion electromagnetic form factor based on 
Particle Data Group p meson parameters, rather than direct data, and they assign a 
very small uncertainty to their evaluation of the dispersion integral in that important 
low-energy region. Comparison of their results with more recent data-driven studies 
suggest that their evaluation of the p contribution was too large. Part of the problem 
may be traced to their use of a relatively low mass for the p, 768.5 MeV, and larger 
width T p ~ 150 MeV when compared with more recent detailed studies of e + e~ 
and tau decay data which suggest m p ~ 776 MeV and V p ~ 146 MeV. At the very 
least, the uncertainty inherent in their approach should have been considerably larger. 
However, if one merely changes the p mass and width to the values currently more 
consistent with data, the AY value of a M (Had; 1) is (very) roughly lowered by ~ 100, 
(xlO -11 ) a significant reduction. To put it on a level with more data-driven recent 
studies, it would need to be fully updated. 

A second more disturbing problem with the AY95 analysis is its use of perturbative 
QCD down to very low energies ~ 1.4 GeV and its unconventional treatment of heavy 
quark thresholds used in conjunction with resonance contributions evaluated by other 
authors. Rather than comment on the details of their analysis, we consider their 
results for another dispersion integral, 

aMl , R(s) 
Aa had = -— ^Re/ ds- -, 14 

the hadronic loop corrections to the fine structure constant Aa^ d (m 2 z ). That quantity 
has been evaluated by many authors fTg ]23|, [Z§ ]3§ ^ |38|, §9] [HJ because 



of its critical use in comparing a, mz and G p with precision measurements of mw, 
sin etc. In Fig.|, we compare various evaluations of that quantity. Notice, its 
significant improvement in the more recent evaluations primarily because of improved 
data. The only calculation which deviates significantly from the others is the value 

A «had( m l) = 289 - 4 ± 4 - 4 x 10 " 4 (AY95) (15) 



given in the text of Adel and Yndurain |yj . (As in the case of hadronic contributions 
to a M , AY95 give very different results for Aaf^ d (m|), called Aa^, in their text and 
abstract.) 

The deviation of Eq.|15] by about 3a from the central value of more up-to-date 
evaluations, which give ~ 276-277 x 10~ 4 would seem to invalidate the perturbative 
analysis of heavy quark thresholds used in conjunction with explicit resonance contri- 
butions for that quantity (suggesting double counting) and thereby call into question 
its use for — 2 studies. We might note that comparison of mw and sin 2 9$ can also 
be used in the Standard Model framework to obtain (a less precise) Aa|^ d (m|) which 
also disagrees with Eq.|i~5|. Again, it seems that one should lower the AY estimate 
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Figure 5: A comparison of the calculated values of Aah a£ z(M§). 



of the high energy hadronic contributions to a M (Had; 1). Taken together with the 
p resonance region reduction mentioned above, such a shift would then move AY's 
central value down to almost exactly DH98b. Of course, a full revision of AY should 
actually be made in any serious update of their approach. 

Based on the above comparisons, we conclude that although CLY85 and AY95 
represent pioneering efforts to incorporate perturbative QCD into a M , their results 
are not state-of-the-art, are not well supported by data and should not be seriously 
considered on a par with more recent evaluations of a M (Had; 1). Furthermore, an 
update of those analyses is likely to significantly reduce their central values to the 
level very similar to more recent evaluations. 

We should also note that several more recent (unpublished) updates or new studies 
of a M (Had; 1) beyond those in Table 1 have appeared. Eidelman and Jegerlehner 
updated their 1995 result with newer data and found [f43[ | 

6967(119) x KT 11 (EJ99) (16) 
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and Jegerlehner has further updated that result (from e + e data) to [44 

6974(105) x KT 11 (J2000) (17) 

A more recent theory-driven analysis (but with parametrization set by e + e~ and r 
decay data) due to Narison found j|5] 



6970(76) x 10~ n (N2001) (18) 

All of these results are in good accord with DH98b which was used in comparison 
with a^ xp (![]. Any of the up-to-date studies will give a deviation between 2 and 2.6 
sigma when compared with experiment. 

5 Summary and Conclusions 

The QED and weak contributions to (<? — 2) are known to an accuracy well below what 
might be accessible to E821. The higher-order hadronic contribution also appears to 
be under control. The recent use of perturbative QCD along with r-decay data have 
greatly improved the uncertainty in the leading hadronic contribution, and super- 
sedes the earlier analyses which only used (more limited) data from electron-positron 
annihilation to hadrons as input. It is clear from the literature that one needs to be 
careful when using parameterizations in the important p region, unless well supported 
by data. The use of perturbative QCD at relatively low energies, which was pioneered 
by CLY85 has led to substantial gains in our understanding, and the agreement of the 



recent BES data [12fl with the QCD predictions of DH98a |23| gives one confidence 
in the validity of the role of modern QCD studies in this discussion. 

The knowledge of the hadronic contribution to has improved dramatically from 
the mid 1980s when the (g — 2) experiment began. New high quality e + e~ data have 
become available from Novosibirsk (y/s = threshold - 1.4 GeV), and from Bejing 
(i/i = 2.0 - 5.0 GeV). The entire sample of LEP r-decay data, as well as the CLEO r 



data [27| have also become available. An updated global analysis of all of these data, 
along with QCD information is underway by Eidelman, Davier and Hocker. These 
authors have now spent many years working on this topic, and we welcome their 
combined efforts to produce a new value for a At (Had; 1). This new more complete 
analysis should make all of the previous analyses obsolete. 

Currently the best published evaluation of a At (Had; 1) is the work by Davier and 
Hocker |53|, [23j. We see no reason to ignore the substantial amount of additional 
information which went into these analyses, and rather to use an earlier evaluation 
which contains only part of our current knowledge. To us, the assertion that one can 
ignore the latest information, or the claim that earlier, and by necessity, less complete 
analyses are on an equal footing is incomprehensible, particularly when the earlier 



approach |2T|] fails so badly for Aahad^l) anc ^ has other serious deficiencies. Both in 
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experimental and theoretical physics, it has been the course of scientific progress to 
have new information and improved analyses of problems replace the results of earlier 
efforts. 

We are fortunate that at this point several independent authors have considerable 
experience in the delicate issues which are involved in obtaining a value for a M (Had; 1). 
They have brought new insights to the field. We anxiously await their updated 
evaluation of this important number. 

We wish to thank Michel Davier for helpful comments on a number of points. LR 
wishes to thank J. Bijnens, D.H. Brown, T, Kinoshita and E. de Raphael for useful 
discussions. This work was supported in part by the U.S. NSF and the U.S. DOE. 
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